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1 Introduction 

The fractional Brownian motion (fBm) B = {Bt, f > 0} of Hurst index H e (0, 1) is a centered 
and continuous Gaussian process verifying Bq ^ a.s. and 

E[iB,-B,r] = \t-s\'" (1) 

for all s,t > 0. Observe that B^^^ is nothing but standard Brownian motion. Equality ([1]) implies 
that the trajectories of B are {H — £)-H61der continuous, for any e > small enough. As the fBm 
is selfsimilar (of index H) and has stationary increments, it is used as a model in many fields 
(for example, in hydrology, economics, financial mathematics, etc.). In particular, the study of 
stochastic differential equations (SDEs) driven by a fBm is important in view of the applications. 
But, before raising the question of existence and/or uniqueness for this type of SDEs, the first 
difficulty is to give a meaning to the integral with respect to a fBm. It is indeed well-known that B 
is not a semimartingale when H 1/2. Thus, the Ito or Stratonovich calculus does not apply to 
this case. There are several ways of building an integral with respect to the fBm and of obtaining 
a change of variables formula. Let us point out some of these contributions: 

1. Regularization or discretization techniques. Since 1993, Russo and Vallois |31| have developed 
a regularization procedure, whose philosophy is similar to the discretization. They introduce 
forward (generalizing Ito), backward, symmetric (generalizing Stratonovich, see Definition [3] 
below) stochastic integrals and a generalized quadratic variation. The regularization, or dis- 
cretization technique, for fBm and related processes have been performed by [T^ [T71 [5^ [55] . 
in the case of zero quadratic variation (corresponding to i? > 1/2). Note also that Young inte- 
grals [3S], which are often used in this case, coincide with the forward integral (but also with 
the backward or symmetric ones, since covariation between integrand and integrator is always 
zero). When the integrator has paths with finite p- variation for p > 2, forward and backward 
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integrals cannot be used. In this case, one can use some symmetric integrals introduced by 
Gradinaru et al. in JT^ (see §2 below). We also refer to Errami and Russo [TT] for the specific 
case where H > 1/3. 

2. Rough paths. An other approach was taken by Lyons [20^ . His absolutely pathwise method based 
on Levy stochastic areas considers integrators having p- variation for any p > 1, provided one 
can construct a canonical geometric rough path associated with the process. We refer to the 
survey article of Lejay [18] for more precise statements related to this theory. Note however 
that the case where the integrator is a fBm with index iJ > 1/4 has been studied by Coutin 
and Qian [7] (see also Feyel and de La Pradelle [13] )• See also Nourdin and Simon [5S] for a 
link between the regularization technique and the rough paths theory. 

3. Malliavin calculus. Since fBm is a Gaussian process, it is natural to use a Skorohod approach. 
Integration with respect to fBm has been attacked by Decreusefond and Ustiinel [8] for H > 1/2 
and it has been intensively studied since (see for instance [1 , 2, 6J ) , even when the integrator is a 
more general Gaussian process. We refer to Nualart's survey article [57] for precise statements 
related to this theory. 

4. Wick products. A new type of integral with zero mean defined using Wick products was intro- 
duced by Duncan, Hu and Pasik-Duncan in [10], assuming H > 1/2. This integral turns out 
to coincide with the divergence operator. In [3], Bender considers the case of arbitrary Hurst 
index H G (0, 1) and proves an Ito formula for generalized functionals of B. 

In the sequel, we will focus - in dimension one - on SDEs of the type: 

r dXt = a{Xt) dBt + b{Xt)dt, t e [0, T] 
\Xo = xoeM 

where cr, & : M — > M are two continuous functions and H £ (0, 1). Our principal motivation is to 
describe a simple theory - from our point of view - allowing to study the SDE (0), for any H g 
(0, 1). It is linked to the regularization technique (see point 1 above). Moreover, we emphasize that 
it is already used and quoted in some research articles (see for example [4l [Ml [2T1 [22l l24l l25l |26] 1 . 
The aim of the current paper is, in particular, to clarify this approach. 

The paper is organized as follows. In the second part, we will consider several definitions 
of solution to ([2]) and for each of them we will study under which condition one has existence 
and/or uniqueness. Finally, in the third part, we will examine whether or not the canonical scheme 
associated to ^ converges, when the integral with respect to fBm is defined using the Russo- Vallois 
symmetric integral. 

2 Basic study of the SDE 

In the sequel, we denote by B a fBm of Hurst parameter H S (0, 1). 

Definition 2.1 Let X,Y he two real continuous processes defined on [0,T]. The symmetric inte- 
gral (in the sense of Russo- Vallois) is defined by 

r Y^d°Xu - hm in prob ^"+" + x ^"+" ~ ^" du, (3) 

provided the limit exists and with the convention that Yt — Yt and Xt — Xt when t > T. 

Remark 2.2 If X, Y are two continuous semimartingales then Yud°Xu coincides with the 
standard Stratonovich integral, see [3T] . 

Let us recall an important result for our study: 

Theorem 2.3 (see /7?| /. p. 793). The symmetric integral f{By)d°Bu exists for any / : R ^ R 

of class if and only if H £ (1/6, 1). In this case, we have, for any antiderivative F of f : 

F{Bt) = F{0)+ [ f[B^)d°Bu. 
Jo 
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When H < 1/6, one can consider the so-called m-order Newton-Cotes functional: 

Definition 2.4 Let / : M" ^ M (with n > 1) be a continuous function, X : [0,T] x Q ^ R 
and Y : [0, T] x i7 — > M" be two continuous processes and m > 1 be an integer. The m-order 
Newton-Cotes functional of (/, Y, X) is defined by 

= liiir in prob / ( / /(F^ + - Yu))i^„r{df3)) 

provided the limit exists and with the convention that Yt — Yt and Xt — Xt when t > T. Here, 
= -\- 5i) and 

Vm^ }. [ I I ■■ 7 «")<5j7(2m-2), m > 2, (4) 

^=0 ^-^0 k^j 
Sa being the Dirac measure at point a. 

Remark 2.5 • The 1-order Newton-Cotes functional f{Yu)d^'-^'^Xu is nothing but the sym- 
metric integral Jq f{Yu)d° X^ defined by On the contrary, when to > 1, the m-order 
Newton-Cotes functional f{Yu)d^^ "'Xu is not a priori a"true" integral. Indeed, its defi- 
nition could be different from jj" /(y„)c?^'-^'™X„ even if f{Y) = f{Y). This is why we call it 
^^functionaF instead of ^HntegraF . 

• The terminology Newton- Cotes functional" is due to the fact that the definition of iy„i via ([4]) 
is related to the Newton-Cotes formula of numerical analysis. Indeed, Vm is the unique discrete 
measure carried by the numbers j /{2m ~ 2) which coincides with Lebesgue measure on all 
polynomials of degree smaller than 2m — 1 . 

We have the following change of variable formula. 

Theorem 2.6 (see \15^ . p. 793). Let m > 1 be an integer. The m-order Newton-Cotes functional 
Iq f{Bu)d^'^'"'Bu exists for any f : R ^ R of class C^^+i if and only if H e (1/(4to + 2), 1). Ln 
this case, we have, for any antiderivative F of f : 

F{Bt) = F(0) + r f[Bu)d'''' "'B^. (5) 
Remark 2.7 An immediate consequence of this result is that 

r f[B^)d'''''^B^ = r /(B„)dNc,n^^ ^ _ ^(0) 

Jo Jo 

when m > n, / is C^™"*"^ and H € (l/(4n -I- 2), 1). Then, for / regular enough, it is possible to 
define the so-called Newton-Cotes functional f{Bu)d^'^Bu without ambiguity by: 

r /(i3„)d^Ci?„ := r if e (l/(4n + 2), 1). (6) 
Jo Jo 

In the sequel, we put uh — inf{n > 1 : H > l/(4n -I- 2)}. An immediate consequence of ([5]) and 

dSl) is that, for any H G (0, 1) and any / : M K of class c4"«+\ we have: 

F{Bt)^F{0)+ r f{B^)d''''B^, 
Jo 

where F is an antiderivative of /. 
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To specify the sense of Jl;a{Xs)dBs in ©, it now seems natural to try and use the Newton-Cotes 
functional. But for the time being we are only able to consider integrands of the form f{B) with 
/ : K ^ R regular enough, see That is why we first choose the following definition for a 
possible solution to 

Definition 2.8 Assume that a e C'*""+^ and that 6 G C°. 

i) Let £i be the class of processes X : [0,T] x ^ R verifying that there exists / ; R — > R 
belonging to C''""^"'^ and such that, for every t G [0, T], Xt = f{Bt) a-s. 
a) A process X : [0, T] x f? — > R is a solution to ^ if: 

• X e€i, 

• yt e [0, T],Xt^xo + Jo aiXs)d''^B, + b{Xs)ds. 

Remark 2.9 Note that the first point of definition ii) allows to ensure that the integral 
J*a{Xs)d^'^Bs makes sense (compare with the adaptedness condition in the Ito context). 

We can now state the following result. 

Theorem 2.10 Let a G C^""+^ be a Lip schitz function, b be a continuous function and xq be a 
real. Then equation 0j admits a solution X in the sense of Definition \2.8\ if and only ifb vanishes 
on S(R), where & is the unique solution to 6' — a o @ with initial value S(0) — xo- In this case, 
X is unique and is given by Xt = &{Bt). 

Remark 2.11 As a consequence of the mean value theorem, 6(R) is an interval. Moreover, it is 
easy to see that either 6 is constant or 6 is strictly monotone, and that inf 6(R) and sup6(R) 
are elements of {ct = 0} U {icx^}. In particular, if a does not vanish, then 6(R) — R and an 
immediate consequence of Theorem 12.101 is that ([2]) admits a solution in the sense of Definition 
[m if and only if & ee 0. 



Proof of Theorem 12.101 Assume that Xt — f{Bt) is a solution to ([2]) in the sense of Definition 
mi Then 

f{Bt) = xo+ [ <Jo f{Bs)d''^Bs + [ bo f{B,)ds = G{Bt) + [ bo f{B,)ds, (7) 
Jo Jo Jo 

where G is the antiderivative oi a o f verifying G(0) — a;o- Set h — f — G and denote by fi* the 
set oi Lo ^ fi such that t i-^ Bt{Lo) is differentiable at least one point t^ G [0,r] (it is well-known 
that P(/?*) = 0). If h'{Bta{uj)) ^ for one (a;, to) € f2 x [0,T] then h is strictly monotone in a 
neighborhood of Btg{uj) and, for \t — to\ sufficiently small, one has Bt{uj) = h~^ {J* b{Xs{Lo))ds) 
and, consequently, ui G f2* . Then, a.s., h'{Bt) = for all t G [0,r], so that ft, = 0. By uniqueness, 
one deduces f = 6. Thus, if ([2]) admits a solution X in the sense of DefInition l2.81 one necessarily 
has Xt = &{Bt). Thanks to 0, one then has b o 6{Bt) = for all t G [0,T] a.s. and then b 
vanishes on S(R). □ 

Consequently, when the SDE Q has no drift 6, there is a natural solution. But what can we do 
when 6^0? 

Denote by jz/ the set of processes A : [Q,T] x f2 ^ R having C^-trajcctories and verifying 
E(e^^o^^='^'*) < oo for at least one A > 1. 

Lemma 2.12 Let A e £/ and to G N*. Then f{B^ + A^)d^'^'"'B^ exists for any f -.R^R of 
class C^™"*"^ if and only if H > 1/ (4m + 2). In this case, for any antiderivative F of f , one has: 

F{Bt + At) = FiAo) + [ /(B„ + A„)d^^'™B„ + / /(B„ + Au)Kdu. 

Jo Jo 
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Proof. Set B ^ B + A. On the one hand, using the Girsanov theorem in 28 and taking into 
account the assumption on A, we have that i? is a fBm of index H under some probabihty Q 
equivalent to the initial probability P. On the other hand, it is easy, by going back to Definition 
\IH to prove that /{B^ + Au)d^'^^"^B^ exists if and only if /{B^ + Au)d^^^"'{B^ + A,,) 
does, and in this case, one has 

r f{B^+A^)d'''''^{B^ + A^)= r f{B^+A^)d'''''^B^+ f{B.^ + A^)A'^du. 
Jo Jo Jo 

Then, since convergence under Q or under P is equivalent, the conclusion of Lcmma l2.12l is a direct 
consequence of Theorem 12.61 □ 



Then, as previously, it is possible to define a functional (still called Newton-Cotes functional) 
verifying, for any H £ (0, 1), for any / : M — > IR. of class 0''""+^ and any process A G 

F{Bt + At) = FiAo) + [ f{B^ + Au)d''^B^ + [ + A„Xdu, 

Jo Jo 

where F is an antiderivative of /. 

Now, we can introduce an other definition of a solution to ([2]): 

Definition 2.13 Assume that a e C^^h+i ^nd that b £ C° . 

i) Let €2 be the class of processes X : [0, T] x J7 — > R such that there exist a function / : R — > R 
in C^"H+i Q^Yid a process A £ such that ^0 = and, for every t G [0,T], Xt = f{Bt + At) a.s. 

ii) A process X : [0, T] x i? ^ R is a solution to ^ if: 

. e [0, TlXt = XQ + a{Xs)d^'^Bs + b{X,)ds. 

Theorem 2.14 Let a G c4nff+i ^ Lipschitz function, b be a continuous function and Xq be a 
real. 

• If ct^xq) ~ then (0j admits a solution X in the sense of Definition \2.13\ if and only if 
b{xo) = 0. In this case, X is unique and is given by Xt = Xq. 

• If a{xf)) 7^ 0, then (0) admits a solution. If moreover infa \cr\ > and b G Lip, this solution is 
unique. 

Proof. Assume that X = f{B + A) is a solution to ^ in the sense of Definition 12.131 Then, we 
have ^ ^ 

f{Bt + At) = G{Bt+At)- j aiXs)A'Js+ [ 6(X,)ds (8) 

Jo Jo 

where G is the antiderivative of cr o / verifying G{0) ~ x^. As in the proof of Theorem 12. 101 we 
obtain that f = & where 6 is defined hy & = a o & with initial value S(0) ~ xq. Thanks to ([8]), 
we deduce that, a.s., we have b o &{Bt + At) — a o &{Bt + At) A't for all t G [0, T]. Consequently: 

• If a{xQ) = then & = xq and b{xo) = 0. 

• If a{xQ) ^ then 6 is strictly monotone and the ordinary integral equation 



At = I ^-—{Bs+ As)ds 



admits a maximal (in fact, global since we know already that A is defined on [0,T]) solution 
by Peano's theorem. If moreover infR \a\ > and b G Lip then — G Lip and A is 
uniquely determined. □ 



The previous theorem is not quite satisfactory because of the prominent role played by xq. That 
is why we will finally introduce a last definition for a solution to ([2]) . We first need an analogue of 
Theorem 12.61 and Lemma [2. 121 
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Theorem 2.15 (see \2^ . Chapter 4)- Let A be a process having -trajectories and m > 1 be 
an integer. If H > l/(2m+ 1) then the m-order Newton- Cotes functional f{Bu, Au)d^'^-^^Bu 
exists for any f : ^ R of class C^™'"'^. In this case, we have, for any function F : —> M. 
verifying F^ = f : 

F{BT,AT)^FiO,Ao)+ [ f{Bu,Au)d'''' "'Bu+ I F^{Bu, Au)Aldu. 

Jo Jo 

Remark 2.16 • F^ (resp. F^) means the derivative of F with respect to a (resp. b). 

• The condition is here H > l/(2m + 1) and not H > l/(4m + 2) as in Theorem 12.61 and 
Lemma Thus, for instance, if A e ^, if 5 : M ^ M is and if /i : ^ M is C^'^ then 

g{Bs + As)d°Bs exists if (and only if) H > 1/6 while h{Bs, As)d°Bs exists a priori only 
when H > 1/3. 

• We define mn — inf{m > 1 : > l/(2m+l)}. As in the Remark l2.71 it is possible to consider, 
for any H G (0,1) and without ambiguity, a functional (still called Newton-Cotes functional) 
which verifies, for any / : — > M of class C^™"^^ and any process A having C^-trajectories: 



F{Bt,At) = F{0, Ao) + / f{Bu, Au)d''''Bu + / F'^{B,,,Au)K 

Jo Jo 



du. 



where F is such that F^' — f. 
Finally, we introduce our last definition for a solution to 
Definition 2.17 Assume that a e C^'"" and b e C°. 

i) Let £3 be the class of processes X : [0, T] x^? — > R verifying that there exist a function / : R'^ — > R 
of class C^™"'^ and a process A : [0,T] x ^? ^ M having -trajectories such that = and 
verifying, for every t £ [0,T], Xt = f{Bt,At) a.s. 
a) A process X : [0, T] x J7 — > R is a solution to if: 

• xe€3, 

• Vi e [0, T],Xt = xo + Jo a(X,)dNC5^ + jt b{X,)ds. 

Tiieorem 2.18 Let a G C^, b be a Lipschitz function and xq be a real. Then the equation |^ 
admits a solution X in the sense of Definition \2.17\ Moreover, if a is analytic, then X is the 
unique solution of the form f(B,A) with f analytic (resp. of class ) in the first (resp. second) 
variable and A a process having -trajectories and verifying Aq — 0. 



Remark 2.19 • If > 1/3, one can improve Theorem l2.18l Indeed, as shown in [26 , uniqueness 
holds without any supplementary condition on a. Moreover, in that reference, another meaning 
to ^ than Definition 12.171 is given, using the concept of Levy area. 

• In 25J, one studies the problem of absolute continuity in equation where the solution is in 
the sense of Definition 1 2. 171 It is proved that, if a{xo) ^ 0, then C{Xt) is absolutely continuous 
with respect to the Lcbesgue measure for all t g]0,T]. More precisely, the Bouleau-Hirsch 
criterion is shown to hold: if xt — xo + Jq b{xs)ds and tx = sup{i G [0,T] : xt ^ IntJ} where 
J = (T~^ {{0}) then C{Xt) is absolutely continuous if and only ii t > t^. 

• We already said that, among the m-order Newton-Cotes functionals, only the first one (that is, 
the symmetric integral, defined by ([3])) is a "true" integral. For this integral, the main results 
contained in this paper are summarized in the following table (where / denotes a regular 
enough function and A a process having C^-trajectories): 
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If we use 
Definition 



[231 



Em 



we have 
to choose 
H e 



(1/6,1) 



(1/6,1) 



(1/3,1) 



X is then 
of the form 



f{B) 



f(B + A) 



f{B,A) 



we have 
existence if 



cr e n Lip, 
G C° and 

&|6(K) = 

(T e n Lip, 
6 e CO + 

i) cr(a;o) = 
b{xo) = 

or 

ii) cr (.Tq) ^ 
and 6 G Lip 



and uniqueness 
if moreover 



*)- 

ii) infffi |(t| > 
and b G Lip 



See 

Theorem 



and [ 



Table 1. Existence and uniqueness in SDE Xt = .tq + / a{Xs)d° Bg + / b{Xs)ds 

Jo Jo 

Proof of Theorem 12.181 Let us remark that the classical Doss-Sussmann [HI [33] method gives a 
natural solution X of the form f{B,A). Then, in the remainder of the proof, we will concentrate 
on the uniqueness. Assume that X — f{B, A) is a solution to ((2|) in the sense of Definition 12.171 
On the one hand, we have 



Xt = Xo 



b{Xs)ds 



f aiX,)d''''Bs+ f 
Jo Jo 

,0+1 <Jof{Bs,As)d''''Bs+ f bof{Bs,As)ds 
Jo Jo 



(9) 



On the other hand, using the change of variables formula, we can write 



Xt^xo 



r/fc(Ss,^)d^^Bs+ f f,{B,,AM'sds. 
Jo Jo 



(10) 



Using © and pI7|) . we deduce that t ^ /J if{Bs, As)d^''^ Bs has C^-trajectories where ip /ft— cro/. 
As in the proof of Theorem 12. 14[ we show that, a.s.. 



Similarly, we can obtain that, a.s.. 



V^G]0,T[, (^(Bt,A) = 0. 



(11) 



VfcGN, V<G]0,T[, ^(B,,^,) = 0. 

If a and /(., y) are analytic, then (^(., y) is analytic and 

G]0, r[, -ix G M, ^(x, At) = j[{x, At)^ao f{x, At) - 0. (12) 
By uniqueness, we deduce 

yt G [0, T], Vx G M, fix, At) = u{x, At), 

where u is the unique solution to uj, = cr(u) with initial value u(0, y) = y for any y G M. In 
particular, we obtain a.s. 



"it e[{),T], Xt ^ f{Bt,At) ^ u{Bt,At). 
Identity ([9]) can then be rewritten as: 



(13) 
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Xt = XQ+ [ aou{Bs,As)d^'^Bs+ [ b o u{Bs, A,)ds, 
Jo Jo 

while the change of variables formula yields: 

Xt=xo+ f u'^{B,,As)dP''Bs+ I u'^{B,,A,)A'Js. 
Jo Jo 



Since u'j^ = a o u, we obtain a.s.: 



Vt e [0,T],bouiBt,At)^u',iBt,At)A't. (14) 
But we have existence and uniqueness in (|14p . Then the proof of Theorem is done. □ 



3 Convergence or not of the canonical approximating schemes 
associated to SDE when d = d° 

Approximating schemes for stochastic differential equations ([2]) have already been studied only in 
few articles. The first work in that direction has been proposed by Lin [TS] in 1995. When H > 1/2, 
he showed that the Euler approximation of equation ([2]) converges uniformly in probability~but 
only in the easier case when a{Xt) is replaced by cr(t), that is, in the additive case. In 2005, I 
introduced in [Mj (see also Talay [M] ) some approximating schemes for the analogue of © where 
B is replaced by a Holder continuous function of order a, for any a € (0, 1). I determined upper 
error bounds and, in particular, my results apply almost surely when the driving Holder continuous 
function is a path of the fBm B, for any Hurst index H G (0, 1). 

Results on lower error bounds are available only since very recently: see Neuenkirch (21) for 
the additive case, and Neuenkirch and Nourdin [22] (see also Gradinaru and Nourdin [M]) for 

equation In [22], it is proved that the Euler scheme X ~ {X associated to ([2]) verifies, 

under classical assumptions on a and b and when i? G (^, 1), that 

^2H-i ^ _x^]^ r a'{X,)DsXids, as n ^ oo, (15) 

2 Jo 

where X is the solution given by Theorem 12.181 and DX its Malliavin derivative with respect 
to B. Still in [22], it is proved that, for the so-called Crank-Nicholson scheme X — {X'^"^}„gN» 
associated to ([2]) with b = and defined by 

^ik+l)/n = ^k/n + U'^i^k/n) + ^(^(fc+l)/™)) ~ -Sfc/„), (16) 

fc e {0, . . . ,n - 1}, 
we have, for a regular enough and when G (|, ^): 

for any a < 3i7 - 1/2, n" ~ Xi] ^ as n ^ oo, (17) 

where X is the solution given by Theorem l2.10l Of course, this result does not give the exact rate 
of convergence but only an upper bound. However, when the diffusion coefficient a verifies 

<t{x)'^ — ax? + /3x + 7 for some a, /?, 7 G M, (18) 

the exact rate of convergence can be derived: indeed, in this case, we have 

^3ff-i/2 _ ^ ^(^^) G, as n ^ 00, (19) 
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with G a centered Gaussian random variable independent of Xi, whose variance depends only 
on H. Note also that, in [TJ, the exact rate of convergence associated to the schemes introduced 
in [21] are computed and results of the type pT|) - P^ are obtained. 

In this section, we are interested in whether scheme (|16p converges, according to the value of 
H and the expression of a. First of all, this problem looks easier than computing the exact rate 
of convergence, as in [HI [22]. But, in these two papers, no optimality is sought in the domain of 
validity of H. For instance, in ([T7|) . we impose that H > 1/i although it seems more natural to 
only assume that H > 1/6. 

Unfortunately, we were able to find the exact barrier of convergence for (|16p only for particular 
a, namely those which verify p8)) . In this case, we prove in Theorem 13.11 below that the barrier 
of convergence is if = 1/6. In the other cases, it is nevertheless possible to prove that the scheme 
(fT6|) converge when iJ > 1/3 (see the proof of Theorem 13. ip . But the exact barrier remains an 
open question. 

The class ^TE\\ is quite restricted. In particular, I must acknowledge that Theorem 13.11 has a 
limited interest. However, its proof is instructive. Moreover it contains a useful formula for ^^"^ 
(see Lemma [3^ . which is the core of all the results concerning the Crank-Nicholson scheme proved 
in [ig (see also [H]). 

Now, we state the main result of this section: 

Theorem 3.1 Assume that a e C^(R) verifies H^) . Then the sequence defined by 

converges in if and only if H > 1/6. In this case, the limit is the unique solution at time 1 to 
the SDE Xt = a{Xg)d° Bg, in the sense of Definition \2.8\ and given by Theorem \2.1(A 

Remark 3.2 When cf{x) = x it is easy to understand why x'^^ converges in \? if and only if 
H > 1/6. Indeed, setting = -B(fe+i)/„ — B^/n, we have 

^1 -^"11 1 _ M« -^oexpO In- 

k=0^ 2^fe 'fc=0 ^ 2^k' 

but 



1 + -A^ ^ 1 

In - + -{AD' + -{Alf + 0{{Alf), 

and, because X]fc=o ^fc ~ -^i ^^d by using Lemma [3.31 below, one has that X^""* converges if and 
only li H > 1/6 and that, in this case, the limit is xq exp(i?i). 

As a preliminary of the proof of Theorem l3.1[ we need two lemmas: 
Lemma 3.3 Let m > 1 be an integer. 

• We have 

n — l ^ 

/ (Bif^iD/j-^ — -Bfc/n)^'" converges m L^ as n ^ oo if and only if H > . 

k=0 

In this case, the limit is zero if H > l/2m and is (2TO)!/(2'"m!) if H = l/2m. 

• We have 

n— 1 ^ 

^^(^(fc+i)/n ~ ^fc/n)^™'''^ converges inlJ^ as n ^ (X if and only if H > -. 

k=Q ^ ^ 

In this case, the limit is zero. 

Proof of Lemma 13.31 The first point is an obvious consequence of the well-known convergence 



n-1 

^2mH-i Y^{B^,+,y„ - Bfc/J^" (2m)!/(2'"m!), as n 



k=0 
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Let us then prove the second point. On the one hand, for H > 1/ (Am + 2), we can prove directly 
that 



E[{Bik+i)/n - - S£/n)""+'] ^ 0, as n ^ oo, 



kl=0 



by using a Gaussian hnear regression, see for instance 16J, Proposition 3.8. On the other hand, it 
is well known that, when H < 1/2, 



n-l 

^(2,n+i)H-i/2 Y^^B^^^iy^ - Bfe/„)""+' ^ N(0, a^,^), as n ^ oo, for some > 

fe=0 

(use, for instance, the main result by Nualart and Peccati [IH]). We can then deduce the non- 
convergence when H < \/{Am + 2) as in 15j, Proof of 2(c), page 796. □ 

Lemma 3.4 Assume that a £ C^(R) is hounded together with its derivatives. Consider (j> the flow 
associated to a, that is, (f>{x, •) is the unique solution to y' ~ cr(y) with initial value y(0) — x. Then 
we have, for any I S {0, 1, . . . , n}.- 

k=0 k=0 k=0 

(20) 

Here we set 

h = ^ — fi = ^ — — and /5 = \ 1 \ \ , 

■' 24 ' 48 80 15 40 20 80 

Zifc = B(fc+i)/„ - Bh/n, when n G N and fc G {0, 1, . . . , n - 1} 

and 

AP{B)^ max \(AlY\ , when p . 

fc— 0,...,n— 1 

Proof of Lemma 13.41 Assume, for an instant, that a does not vanish. In this case, (f){x, ■) is a 
bijection from M to himself for any x and we can consider 1^9(2;, •) such that 

Vx, t e R : Lp{x, cj)(x, t)) = t and (/)(x, ip{x, t)) = t. (21) 

On the one hand, thanks to (PTjl . it is a little long but easy to compute that 

ip{x,x) =0, 

^'t{x,x) ^l/a{x), 

^'l,{x,x) ^[^a'/a^]{x), 

^f2{^,x) =[(2a'2-aa")M(x), 

^it\(x, x) ^ [(-6a'3 + &aa'a" a^a'")la^] [x) 

^ftlttix, x) = [(24a'4 - 36aa'2c7" + Sa^aV'" + Ga^a''^ - a^^^) / a%x) . 
Then, for u sufficiently small, we have 

<^{x,x + u) - Ux)u - ^{x)u^ + ^"'\-r" {xW + (x)u^ 



2^V-^^" -f- ^ 1.^/" -r 24^ 

" - - -i I A\ 

(x)m5 + 0(1*6). 



On the other hand, using (I16p and some basic Taylor expansions, one has for fc G {0, 1, n — 1}: 

^(fc+l)/« - ^kjn^ ^y^kjJ'^k + — i^fc/„n^fe) H 4 l^fe/nA^fc) 

+ + 8 + —vr)^^kln)\Ak) 

+ \— + 8 + 6 + + -^)y^kln)\^k) 

+ 0{A''iB)). 
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Then, we have 

= A- + UxffjiA-)^ + UixffjiA-r + Mxi-l){A-r + 0{A%B)). 
We deduce, using ([2T|) : 

= + M^ifji^k? + .U{xl%)iA-,r + h{x]^J{Alf + 0{A\B))) . 

Finally, by using the semi-group property verified by 0, namely 

Vx, s, t e M : (f>{<j)ix, t), s) = (j){x, t + s). 

we easily deduce (1^ . 

In fact, we assumed that a does not vanish only for having the possibility to introduce ip. But 
(1^0)1 is an algebraic formula then it is also valid for general a, as soon as it is bounded together 
with its derivatives. □ 

Proof of Theorem 13.11 Assume that a verifies . Although a is not bounded in general, it 
is easy to verify that we still have 0{nA^{B)) as remainder in ([^0]) . Moreover, simple but tedious 
computations show that we can simplify in (|20p to obtain 

= 0(.To, Si + % EVoi^kf + m EVoi^kf + 0{nA%B))). 

Thus, as a conclusion of Lemma 13.31 we obtain easily that X^^ converges to ^(a;o, Bi) if and only 
ifi7>l/6. □ 

Acknowledgement. I am indebted to the anonymous referee for the careful reading of the original 
manuscript and for a number of suggestions. 
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